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Abstract  ^  I 

In  this  paper,  we  present  new  results  on  the  reconstruction  of  signals  from  only 
the  sign  of  the  real  part  of  the  Fourier  transform.  Specifically,  w«  develop  new 
theoretical  results  which  state  conditions  under  which  two-dimensional  signals  are 
uniquely  specified  to  within  a  scale  factor  with  this  information  and  show  that  tjiese 
conditions  include  a  broad  class  of  signals.  Furthermore.  we  apply  this  result,  Yo  the 
problem  of  reconstructing  two-dimensional  signals  from  their  zero  crossings.  Wo  also 
>  present  two  algorithms  for  reconstructing  a  signal  from  sign  information  in  either  the 
time  or  frequency  domain.  ' ..  _ 


L  Introduction 


Signal  reconstruction  from  partial  Fourier  domain  information  has  been  of 
interest  to  a  number  of  different  authors  both  for  particular  applications  and  for  its 
inherent  theoretical  value  [1].  Previous  work  in  this  area  has  involved  developing  con¬ 
ditions  under  which  signals  are  uniquely  specified  with  Fourier  transform  magnitude  or 
phase  [23,4]  or  signed-magnitude  [5]  information  and  developing  practical  algorithms 
for  recovering  signals  from  this  information.  In  this  paper,  we  consider  the  problem  of 
reconstructing  signals  from  only  Fourier  transform  (or  inverse  Fourier  transform)  sign 
information.  This  sign  information  can  be  viewed  as  one  bit  of  the  Fourier  transform 
phase,  without  any  magnitude  information.  Alternatively,  this  information  can  be 
viewed  as  the  zero  crossings  of  the  real  part  of  the  Fourier  or  inverse  Fourier 
transform.  This  latter  viewpoint  suggests  a  number  of  practical  applications  pertaining 

"This  work  has  been  supported  in  part  by  the  Advanced  Research  Projects  Agency  monitored  by 
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to  the  reconstruction  of  signals  from  zero  crossing  information.  One  such  application 
occurs  when  a  signal  is  clipped  or  otherwise  distorted  in  such  a  way  as  to  preserve  the 
zero  crossing  or  level  crossing  information,  and  it  is  desired  to  recover  the  original  sig¬ 
nal.  Another  application  occurs  in  the  theory  of  vision  where  studies  have  stressed  the 
importance  of  edge  detection  as  a  means  of  classifying  and  identifying  images  but  have 
not  succeeded  in  developing  a  strong  theoretical  basis  for  this  work  [6].  A  third  appli¬ 
cation  occurs  in  some  design  problems  where  one  might  want  to  specify  a  filter 
response  (7]  or  antenna  pattern  [8]  in  terms  of  zero  crossing  or  null  points  (such  as  for 
interpolation)  and  derive  the  remainder  of  the  response  from  these. 

In  our  previous  work,  we  have  established  the  importance  of  the  sign  of  the  real 
part  of  the  Fourier  transform  (also  referred  to  as  one  bit  of  Fourier  transform  phase) 
both  experimentally  and  theoretically  [1].  On  the  experimental  side,  we  have  shown 
with  images  that  if  the  correct  one  bit  of  phase  is  combined  with  a  unity  or  average 
magnitude,  the  resulting  image  maintains  many  of  the  features  of  the  original  image, 
and  in  fact,  is  identical  to  the  phase-only  version  of  the  even  (symmetric)  component 
of  the  image.  Furthermore,  if  an  image  is  synthesized  from  the  correct  magnitude  and 
phase  with  the  most  significant  bit  of  the  phase  randomized,  the  result  is  unintelligible. 
On  the  theoretical  side,  we  have  demonstrated  that  much  stronger  results  can  be 
obtained  for  unique  specification  of  signals  with  signed  Fourier  transform  magnitude 
than  with  Fourier  transform  magnitude  alone. 

We  have  also  previously  established  conditions  under  which  a  one-dimensional 
signal  is  uniquely  specified  to  within  a  scale  factor  with  Fourier  transform  sign  infor¬ 
mation  [1],  although  these  conditions  are  fairly  restrictive.  One-dimensional  signals  are 
uniquely  specified  with  this  information  only  if  the  real  part  of  the  Fourier  transform 


contains  a  sufficient  number  of  sign  changes,  or  zero  crossings.  In  the  two-dimensional 
case,  a  'zero-crossing*  is  actually  a  contour  in  the  uj)  plane  and  thus  consists  of  an 
infinite  number  of  points.  Thus,  it  is  reasonable  to  expect  that  in  the  two-dimensional 
problem,  sign  information  in  the  Fourier  domain  is  more  likely  to  uniquely  specify  a 
signal  than  in  the  corresponding  one-dimensional  problem. 

In  this  paper,  we  present  new  results*  on  the  unique  specification  of  two- 
dimensional  signals  with  Fourier  transform  sign  information.  Specifically,  we  develop 
conditions  under  which  a  two-dimensional  signal  is  uniquely  specified  to  within  a  scale 
factor  with  sign  information  alone.  Furthermore,  we  show  that  these  conditions 
include  a  broad  class  of  signals.  These  results  are  presented  in  section  2.  We  also 
develop  a  number  of  extensions  to  these  results  and  discuss  the  problem  of  unique 
specification  with  sign  information  available  only  at  a  finite  set  of  discrete  frequencies. 

In  section  3,  we  apply  our  results  to  the  problem  of  recovering  a  two-dimensional 
signal  from  its  zero  crossings  and  more  generally,  from  its  threshold  crossings.  A  con¬ 
siderable  amount  of  research  has  been  devoted  to  the  problem  of  reconstructing  one¬ 
dimensional  signals  from  zero  crossings  [1344]  and  much  less  to  the  corresponding  two- 
dimensional  problem.  The  problem  of  reconstruction  from  zero  crossings  is  a  dual  to 
the  problem  of  reconstruction  from  one  bit  of  phase  since  it  involves  recovering  a  sig¬ 
nal  from  sign  information  in  the  signal  domain  rather  than  in  the  Fourier  domain  and 
thus  the  results  developed  in  section  2  are  directly  applicable.  In  section  4,  we  discuss 
two  potential  algorithms  for  recovering  a  signal  from  sign  information  in  either 
domain. 


*Sotne  of  the  results  presented  here  were  first  presented  at  the  Inti  Conf.  on  Acoustics,  Speech,  and 
Signal  Processing  in  March  1964  [9], 
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2.  Unique  Specification  with  Fourier  Tramform  Sign  Information 


In  this  section  we  present  new  theoretical  results  on  the  unique  specification  of 
two-dimensional  signals  with  Fourier  transform  sign  information.  After  introducing 
some  notation,  we  will  present  a  number  of  results  which  apply  when  the  sign  informa¬ 
tion  is  available  at  all  frequencies  and  then  extend  these  results  to  situations  where  the 
sign  information  is  available  only  at  a  discrete  set  of  frequencies. 

Notationally,  we  will  use  x  [nj,  nj  to  denote  a  two-dimensional  discrete-time 
sequence,  X  (zj,  zj)  to  denote  its  z-transform,  and  X  (aij,  u>j)  to  denote  its  Fourier 
transform,  ir  - 


*i  "a 

x  <«..  -2)  -  2  2  *  »d  «*•%■*■* 

•* 

The  Fourier  transform  sign  information,  or  one  bit  of  phase,  will  be  defined  as: 


|  1  if  Re  fa,  wj)}  2  0 

otherwise  () 

We  will  also  refer  to  the  even  (symmetric)  and  odd  (antisymmetric)  components  of  a 

signal,  defined  as: 


,  *[|»1,|»2]  +x[-W1,-«»2] 

X0lnUn2l  31  5 


(3) 


Similarly,  Xt(zitz  j)  will  denote  the  z-transform  of  xt  [/i  lt  n  J.  We  will  refer  to  z- 
transforms  as  symmetric  if  they  correspond  to  symmetric  sequences,  that  is,  X  (z  j,  z  2)  »* 
symmetric  if  X  (zlfz 2)  (z^.zj*1).  A  factor  of  X  (zlfz 2)  will  be  said  to  be  a  real 

symmetric  factor  if  it  is  symmetric  as  defined  above  and  if  all  of  its  coefficients  are 


% 


•  5  • 


real.  This  does  not  imply  that  X  (zj,z  2)  will  take  on  only  real  values  but  does  imply 
that  X  (2 !, r  j)  will  be  real  on  the  unit  surface  (zjl  ■  |z2|  ■  L  Furthermore,  the  set  of 
real  symmetric  factors  of  a  z-transform  includes  all  possible  factors  which  satisfy  the 
definition  above  and  is  not  limited  to  irreducible  factors. 


21  Unique  Specification  with  S,^,  wj) 

In  this  section,  we  will  develop  results  on  the  unique  specification  of  signals  with 
•  the  Fourier  transform  sign  information  Sg  («lf  <02)1  when  this  information  is  available 

for  all  frequencies.  These  results  can  be  stated  in  a  number  of  different  forms  since  it 
is  possible  to  uniquely  specify  a  signal  with  Fourier  transform  sign  information  under  a 
number  of  different  sets  of  constraints.  After  stating  and  proving  our  primary  result, 
we  develop  a  number  of  extensions  for  different  types  of  sequences  and  different 
£  definitions  of  Sx  (cu1,  We  also  show  that  the  results  developed  here  apply  to  a 

broad  class  of  two-dimensional  signals. 


21L  Primary  Result 

Since  our  uniqueness  theorem  relies  primarily  on  a  well-established  result  from 
algebraic  geometry,  we  shall  first  state  this  result  without  proof.  The  detailed  proof  is 
available  in  references  [10]  and  [11]. 


Theorem  1  [10,11].  If  X(zltz£  and  Y(zi,z£  are  two-dimensional  poly¬ 
nomials  of  degrers  r  and  x  with  no  common  factors  of  degree  >  0,  then 
there  are  at  most  rx  solutions  to  the  following  equations: 

X(*i.*t)«0  (4) 


In  this  theorem,  the  degree  of  a  polynomial  in  two  variables  is  defined  in  terms  of 
the  sum  of  the  degrees  in  each  variable  (for  each  term),  that  is,  the  degree  of  a  two- 
dimensional  polynomial  p(x,y)  is  equivalent  to  the  degree  of  the  one-dimensional 
polynomial  p(x  ,x).  For  example,  a  two-dimensional  sequence  with  support  over 
0  s  n  lt  ij2  would  have  a  z-transform  which  is  a  polynomial  of  degree  2N . 

Essentially,  Theorem  1  places  an  upper  bound  on  the  number  of  points  where  two 
two-dimensional  polynomials  can  both  be  zero  if  they  do  not  have  a  common  factor. 
As  we  state  and  prove  in  Theorem  2,  this  result  can  be  applied  directly  to  the  problem 
of  unique  specification  of  signals  with  Sz  (ialt  wj).  Consider,  for  example,  two  signals 
x  [" i*  "d  and  y  [R i»  "J  f°r  which  Sx  (<i»lf  wj)  =  Sy  (uj,  First  of  all,  we  note  that  since 
the  real  part  of  the  Fourier  transform  only  contains  information  about  the  even  com¬ 
ponent  of  the  sequence,  we  must  require  that  x  [jij,  nj  be  even  or  be  defined  only  over 
a  nonsymmetric  half-plane*  so  it  can  be  recovered  from  its  even  part.  Also,  we  note 
that  if  Re  (X  (<uj,  wj)}  >0  for  all  (wj,  wj),  then  we  could  not  expect 
sign  (Re  {X  (o>1, 012)})  to  be  sufficient  to  reconstruct  the  original  signal.  Thus,  we  will 
also  assume  that  Re  {X  (o>j,  ojj)}  and  Re  {X  («j,  0*2)}  are  positive  in  some  regions  of  the 
(<■»!,  <1*2)  plane  and  negative  in  other  regions.  The  boundary  between  these  regions  is  a 
contour  where  Re  {X  (u>j,  wj)}  =*  Re  {X  (wj,  <02)}  **  0.  With  some  algebra,  we  can  show 
that  this  implies  two  polynomials  in  (*1,22)  ue  zero  ovcr  the  same  infinite  set  of 
points,  and  thus  from  Theorem  1,  these  polynomials  must  contain  a  common  factor.  If 
we  also  assume  that  these  polynomials  are  nonfactorable,  then  they  must  be  equal  to 
within  a  scale  factor,  and  x  [n  x,  a 2]  and  y  [« 1*  must  be  equal  to  within  a  scale  factor. 

*A  nonsymmetric  half-plane  (NSHP)  region  of  support  is  defined  to  mean  that  if  («,,  *j)  is  in  the  re¬ 
gion  of  support,  then  -*ij)  is  not  in  the  region  of  support  unless  n,  -  n,  -  0. 


Specifically,  we  state  the  following  theorem: 

Theorem  2.  Let  xl/ii./ij  and  y  [« i, « 2]  be  real  two-dimensional 
sequences  with  support  over  a  finite  non-symmetric  half-plane,  with 
Sx  («,,  wj)  *  Sy  («,.  »j).  If  Re  (X  («!,  «2)}  takes  on  both  positive  and  nega¬ 
tive  values  and  Xt(zuz^)  and  Yt(zuz^  are  nonfactorable,  then 
x  [n  1,  *  cy  {i»  1,  nj  f°r  some  positive  constant  e . 

Proof:  From  the  theorem  statement,  we  know  that  Sx  (wj,  =  Sy  («x,  0*2);  we  will 
show  that  xlnj.nJ  =cy  for  some  c.  Since  we  know  that  Re  {X  takes 

on  positive  and  negative  values,  there  must  be  some  region  of  the  (ut,  0*2)  plane  where 
Re  {X  <*>2)}  >  0  and  another  region  where  Re  {X  (<olf  **2)}  <  0.  The  boundary 
between  these  regions  is  a  contour  where  Re  (X  («j,  <02)}  =  0,  or  equivalently, 
X#  («j,  «2)  *0.  Thus, 

*.(*..^1,  ,^,^  -  «  (3) 

if  <dj  and  «2  are  on  this  contour.  Since  Sx  (wj,  <1*2)  *  Sy  («j,  0*2)  for  all  («lf  M2),  equation 
(5)  alsc  holds  for  Yt(zuz£.  Thus,  we  have  an  infinite  set  of  points  where 

X,(z  1**2)  88,0  W 

Since  x,[bi,  bJ  is  nonzero  for  positive  and  negative  values  of  »j  and  b2,X,(zj,  12)  is  a 

polynomial  in  the  variables  Zj,  z2,  Zj_l,  and  *2-1-  However,  if  x  [bj,  bJ  and  y  [b^  bJ 

N  jy 

have  finite  support,  then  we  can  find  integers  N 1  and  iV 3  sucb  that  zj  1  z2  1  X#(zj, zj) 
and  m?'  z2**  (z  1, z 2)  are  polynomials  in  only  Z!  and  z2.  Furthermore, 

,, ’S/’ *.(»„*,) -0  P) 

over  the  contour  in  the  (u>i,  plane  where  Re  {X  (wj,  <*>2)}  -  0,  and  therefore,  over  an 
infinite  set  of  points.  Thus,  by  Theorem  ltX,(zltzj)  and  Yt(zuz2)  must  have  a  com- 


mon  factor.  If  furthermore,  we  assume  that  Xt  (z1,z2)  and  Yt(zuz^)  are  nonfactorable, 
then  Xt(zx,z£  =  clr#(z|fZ2),.x,lt*i»»2]  =cy»l«i»«2]»an^Jf[,,i*,,2l  =cy  [nlf  bJ- 

Instead  of  constraining  the  signals  to  be  real  and  to  have  support  over  a  nonsym- 
metric  half -plane,  we  could  just  as  easily  permit  complex  signals  or  signals  symmetric 
about  the  origin.  We  will  state  one  alternate  form  of  Theorem  2  as  a  corollary  since  it 
will  be  needed  later  in  this  paper  (see  Appendix  1  for  proof): 

Corollary  L  Let  x  [/lj,  nj  and  y[nun^\  be  complex  two-dimensional 
conjugate-symmetric*  sequences  with  finite  support  with 
Sx  (<■>!,  <i>2 )  mSy  (a)},  <i>2).  If  X  («!,  <1)2)  takes  on  both  positive  and  negative 
values  and  X  (zlfz 2)  and  Y(zitz£  are  nonfactorable,  then 
x[nun^\-cy  [nt, nj  for  some  positive  constant  c . 

Note  that  in  this  case  X  («t,  0*2)  is  real  and  that  Corollary  1  directly  constrains 

X  (u>j,  0*2)  and  X  (zlfz  2)  rather  than  the  corresponding  transforms  of  the  even  com- 

« 

ponents. 

2JL2.  Extensions 

Although  Theorem  2  states  a  particular  set  of  constraints  under  which  a  signal  is 
uniquely  specified  with  Sx  (uj,  <■> 2),  a  number  of  different  sets  of  constraints  are  possi¬ 
ble.  In  this  section,  we  shall  introduce  three  extensions  to  Theorem  2.  The  first  result 
imposes  most  of  its  constraints  on  only  one  sequence,  as  opposed  to  imposing  a  number 
of  constraints  on  both  sequences  x[nitn  2]  and  y  [«j,  n^}-  This  result  will  be  convenient 
to  use  when  discussing  algorithms  for  reconstruction  from  Sx  («j1(  <1*2).  The  second 
result  generalizes  the  concept  of  Fourier  transform  sign  information  to  include  a 

*A  sequence  x  [n,,  nj  is  conjugate  symmetric  if  *[<•,,  nj  **•[-*,,  — wj. 
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broader  class  of  definitions  than  that  given  in  equation  (2).  The  third  result  extends 
Theorem  2  to  include  sequences  with  factorable  z-transforms  by  replacing  the  nonfac- 
torability  constraint  with  a  constraint  on  each  factor. 

Let  us  start  by  considering  a  case  where  it  would  be  convenient  to  have  a  result 
which  imposes  most  of  its  constraints  on  only  one  sequence.  Suppose  a  sequence 
x  [/ij,  nj  is  known  to  satisfy  the  constraints  of  Theorem  2,  and  we  would  like  to 
develop  an  algorithm  to  recover  x  [n  j,  nj  from  Sx  (u1(  uij).  It  would  be  convenient  to 
have  a  set  of  constraints  which  guarantee  that  there  are  no  other  sequences  y  [nlt  nj 
with  Sx  (<*>!,  a^)  =  Sy  («!>!,  (i>i)  whether  or  not  y  [/ij,  nj  satisfies  the  constraints  of 
Theorem  2.  While  it  is  simple  enough  to  guarantee  that  a  recovered  sequence  y 
has  the  correct  region  of  support,  it  is  extremely  difficult  to  guarantee  that  Yt(zltz£  is 
nonfactorable.  Without  this  restriction,  l^Czj.zj)  might  contain  a  real  symmetric  fac- 
tor  which  is  positive  on  the  unit  surface,  and  thus  we  could  have 
5,  (wj,  «2)  =  Sy  (<Oj,  wj)  but  x  [«lf  n  J  *  cy  [nlf  nj-  To  avoid  this  problem,  we  note  that 
since  multiplication  by  an  additional  factor  in  the  z -domain  corresponds  to  a  convolu¬ 
tion  in  the  space  domain,  the  sequence  y  [/i  j,  n  J  would  then  have  a  larger  region  of 
support  than  x  [jif,  nj-  Thus,  if  the  exact  size  of  the  region  of  support  of  x  [/ilf  nj  is 
known,  then  this  information,  together  with  Sx  (wlf  taj),  is  sufficient  to  uniquely  specify 
x  [/ij,  nj-  Specifically  (see  Appendix  1  for  proof): 

Theorem  3.  Let  xln^nj  aud  y[«i,  nj  be  real  two-dimensional 
sequences  with  support  over  a  finite  non-symmetric  half-plane  contained  in 
-N  Snltn2^N  with  Sx  (<dj,  torf  *  $7  (iat,  <i>2).  If  x[Y,/V]#0, 

Re  {X  (<■)},  <*>2)}  takes  on  both  positive  and  negative  values,  and  X((z1,z2)  ■* 
nonfactorable,  then  x  [n  lt  nj  =  cy  [a  1,  "J  for  some  positive  constant  c . 


'JVif.Y.d 


Thus,  if  a  sequence  x  [nlf  nj  satisfies  the  constraints  of  Theorem  3,  then  we  know 
that  there  are  no  other  sequences  yliS},^  with  the  same  region  of  support  and  the 
same  Fourier  transform  sign  function.  Therefore,  in  a  somewhat  more  general  sense 
than  was  possible  with  Theorem  2,  we  can  say  that  a  sequence  satisfying  the  con* 
attaints  of  Theorem  3  is  uniquely  specified  to  within  a  scale  factor  with  the  Fourier 
transform  sign  information  and  the  known  region  of  support. 

It  is  also  possible  to  generalize  Theorems  2  and  3  to  allow  unique  specification 
with  broader  classes  of  sign  information  than  the  Sg  uj)  as  defined  in  equation  (2). 
Since  Sz  («,,  wj)  can  be  viewed  as  one  bit  of  Fourier  transform  phase,  we  could  gen¬ 
eralize  Sx  (<«>!,  <1)2)  to  allow  quantizing  the  phase  in  different  ways.  Specifically,  we 
could  define: 

* 

1  if  o  -  s  s  a  +  -• 

-1  otherwise 

or  equivalently, 

s/fo,  «2)  =  sign  (Re  {X  (<ou  <*2)  eJa})  (9) 

The  case  a  -  0  corresponds  to  the  definition  of  Sg  (uj,  uj)  given  in  equation  (2).  Alter¬ 
natively,  since  5,  (mj,  <02)  can  be  viewed  as  the  zero  crossings  of  the  real  part  of  the 
Fourier  transform,  we  could  generalize  Sx  (uj,  wj)  to  allow  crossings  of  an  arbitrary 
threshold  as  follows: 

Sf («i,  Wj)  =  sign  (Re  {X  (uj,  »2)}-0)  (10) 

To  develop  a  result  on  unique  specification  with  generalized  sign  information,  we 
will  combine  these  two  ideas  and  define: 


Sg*  (»i,  <*>2)  *  sign(Re  {If  («,,  -  0) 


(11) 


We  can  then  develop  a  result  similar  to  Theorem  2  for  this  definition  of  sign  informa¬ 
tion  (see  Appendix  1  for  proof): 


Theorem  4.  Let  x[*i,  *2]  and  y  [n  j,  nj  be  real  two-dimensional 
sequences  with  support  over  a  finite  non-symmetric  half-plane,  with 
(«*i» “2)  * «2)  *or  any  a  and  P  such  that  Re {X (ult - p 
takes  on  both  positive  and  negative  values.  Also,  let: 


,  x[nltn^eJa  +x*[-ni,-n^e~Ja  _ 

x[ni,«2J  - - 2 - P  8[ni,  nj 


(12) 


y[»i.«2]  = 


y 


-pBtnj./iJ 


where 


1 1  if  (i*lti»2)  8  (0,0) 

8(«  1,1*2]  =|o  otherwise 

If  X(zltZ2)  and  Yiz^z?)  are  nonfactorable,  then  x [«j, nj  =cy[«1,n2]  for 
(1*1,  "2)  *  (0, 0),  and  x  [0, 0]  cos  a  —  P  ■  e(y  [0, 0]  cos  a  -  P)  for  some  posi¬ 
tive  constant  c . 


The  ambiguity  at  (1*1,112)  8  (0,0)  is  not  just  a  scale  factor;  it  is  a  scaling  of  the 

value  with  respect  to  a  threshold  — Note  that  if  a  =  for  some  odd  integer  k , 

cosa  Z 

then  5xa*  (<■>},  gjj)  contains  no  information  about  x  [0, 0],  and  even  if  p  =  0,  x  [0, 0]  can¬ 
not  be  recovered. 

A  further  extension  to  Theorem  2  involves  replacing  the  nonfactorability  con¬ 
straint  with  a  constraint  on  each  factor.  Let  us  express  Xt(zuz 2)  as  a  product  of  real 
symmetric  factors  F^z^zj).  Observe  that  if  F{(zitzj)=0  for  any  i,  then 
Xt(z  1,  Z2)  =0;  similarly,  if  ^(zj.zj)  =0,  then  at  least  one  of  the  factors  Ft(z1,z2) 
must  be  zero.  Thus,  if  each  factor  contributes  a  set  of  zero  crossing  contours,  each  fac¬ 
tor  will  be  uniquely  specified  by  its  own  zero  crossing  contours,  and  thus  we  can 
develop  a  set  of  conditions  under  which  Xt(zuz£  will  be  uniquely  specified  by  the 


complete  set  of  zero  crossing  contours.  Specifically,  we  state  (see  Appendix  1  for 
proof): 

Theorem  5.  Let  x[ji1,n2]  and  y  [n lf  n J  be  real  two-dimensional 
sequences  with  support  over  a  finite  non-symmetric  half -plane,  with 
S,  («i»  “2)  ■  sy  (“l*  "2)*  Consider  the  factorization  of  and 

Y t  (2  j,  22)  into  real  symmetric  factors  which  are  irreducible  over  the  set  of 
real  symmetric  factors.  If  each  of  these  factors  has  multiplicity  one  and 
takes  on  both  positive  and  negative  values  on  the  unit  surface,  then 
x  l"  i»  "d  *  <7  I™  1*  *2]  for  some  positive  constant  e . 

Theorem  5  can  be  easily  modified  to  permit  complex  signals  or  symmetric  signals, 
to  impose  most  of  its  constraints  on  one  sequence,  or  to  permit  different  definitions  of 
the  Fourier  transform  sign  information  Sx  (wj,  0^)  as  was  done  earlier  in  developing 
extensions  to  Theorem  2.  This  theorem  could  also  be  stated  in  a  slightly  different 
manner  by  considering  all  possible  factorizations  of  X9(zlfz£  and  T, (z  1, * 2)  rather 

1 

than  one  particular  factorization.  In  this  case,  the  requirement  would  be  that  every 
possible  factor  of  X,(z\,z£  and  Yt(zlfz£  must  take  on  both  positive  and  negative 
values  on  the  unit  surface.  The  multiplicity  constraint  is  then  unnecessary  since  if  a 
factor  F  (z  i,z  2)  occurs  with  multiplicity  two  (or  higher),  then  there  will  also  be  a  fac¬ 
tor  F2(zltzj)  which  is  nonnegative  on  the  unit  surface  and  violates  the  constraints  of 
the  theorem. 

Theorem  5  is  related  to  a  result  developed  in  [1]  on  the  unique  specification  of 
one-dimensional  signals  with  Fourier  transform  sign  information.  The  one-dimensional 
result  states  constraints  on  the  zeros  of  the  z-transform  X(z)  which  guarantee  that 
X,(z )  has  zeros  only  on  the  unit  circle  (and  that  these  zeros  are  simple  zeros)  and  thus 
that  sign  {Xe  («)}  is  sufficient  to  uniquely  specify  x[n].  This  result  could  be  stated  in  a 
form  similar  to  Theorem  5  by  considering  the  factorization  of  X,(z )  into  real 


symmetric  factors  irreducible  over  the  set  of  real  symmetric  factors.  In  this  case,  the 
factors  tv  ill  all  be  second-order  and  will  have  either  two  complex  conjugate  zeros  on 
the  unit  circle  or  two  real  zeros.  Restricting  each  of  these  factors  to  have  simple  zeros 
on  the  unit  circle  is  then  equivalent  to  requiring  that  each  factor  take  on  both  positive 
and  negative  values  on  the  unit  circle,  as  in  Theorem  5.  In  addition,  we  require  each 
factor  to  have  multiplicity  one  so  that  we  can  guarantee  that  there  is  a  sign  change  in 
Xt(to)  corresponding  to  each  zero  on  the  unit  circle.  The  primary  difference,  then, 
between  the  one-dimensional  problem  and  the  two-dimensional  problem  is  not  in  the 
type  of  constraints  imposed  but  in  the  likelihood  of  a  signal  satisfying  these  con¬ 
straints.  While  only  a  small  class  of  one-dimensional  signals  will  satisfy  the  appropri¬ 
ate  constraints,  we  will  show  in  the  next  section  that  a  broad  class  of  two-dimensional 
signals  will  satisfy  these  constraints. 

« 

2JL3.  Applicability 

Having  established  a  set  of  conditions  which  guarantee  that  a  signal  is  uniquely 
specified  by  some  partial  information,  it  is  worthwhile  to  determine  whether  or  not 
these  conditions  are  likely  to  apply  to  a  typical  sequence  encountered  in  practice. 
While  it  is  difficult  to  answer  this  question  without  making  some  assumptions  about 
the  type  of  application  involved,  it  is  still  possible  to  show  that  a  broad  class  of  signals 
will  satisfy  the  constraints  of  the  results  developed  earlier.  In  this  discussion,  we  will 
refer  primarily  to  Theorem  3  since  this  is  the  basic  result  which  states  constraints  pri¬ 
marily  upon  the  actual  signal  we  are  trying  to  recover;  similar  results  can  be  easily 
developed  for  the  other  extensions  developed  in  the  preceding  section. 


Let  us  first  note  that  in  many  applications,  the  assumed  region  of  support  of  a 
finite-length  signal  is  somewhat  arbitrary  and  thus  signals  can  often  be  considered  to 
have  any  desired  region  of  support.  However,  since  the  real  part  of  the  Fourier 
transform  is  modified  when  the  assumed  region  of  support  of  a  sequence  is  modified, 
care  must  be  taken  to  ensure  that  the  data  available  as  Sx  (wj,  uj)  truly  corresponds  to 
a  sequence  with  a  region  of  support  as  specified  in  Theorem  3.  In  the  remainder  of 
this  discussion,  we  will  assume  that  the  region  of  support  constraint  is  satisfied. 

With  some  informal  arguments,  we  can  next  show  that  the  probability  of  a  ran¬ 
dom  signal  satisfying  the  remaining  constraints  of  Theorem  3  very  rapidly  approaches 
unity  as  the  number  of  points  in  the  signal  increases.  First,  we  note  that  since  'almost 
all”  two-dimensional  polynomials  are  nonfactorable,  the  nonfactorability  constraint  is 
satisfied  with  probability  one  [442].  Similarly,  almost  all  sequences  will  have 
x  [/V  ,N  ]  *  0.  Next,  consider  the  condition  requiring  Re  [X  (uv  u^)}  to  change  sign.  If 
the  coefficients  of  the  sequence  were  random,  there  would  be  some  small,  but  finite, 
chance  that  the  first  coefficient,  x  [OjO],  would  be  greater  than  the  sum  of  the  magni¬ 
tudes  of  the  others,  and  thus  Sx  («et,  <02)  ”  sign  (x  [OjOD  for  all  (wj,  <*2)-  Thus,  we  cannot 
claim  that  almost  all  sequences  satisfy  the  sign-change  constraint,  but  we  can  argue 
that  the  probability  of  a  random  first-quadrant  sequence  satisfying  this  constraint  very 
rapidly  approaches  unity  as  the  number  of  points  in  the  signal  increases. 

To  see  this,  assume  that  x[n1,/i2]  is  a  first  quadrant  sequence  with  support  N  jXW 2 
and  form  the  one-dimensional  sequence 

a [n 2  +N1/I2]  (13) 

If  A  (2)  has  at  least  one  zero  outside  the  unit  circle,  then  Re  {A («)}  is  guaranteed  to 

have  at  least  one  sign  change  over  the  interval  (0,ir)  [1].  Since  A  (u>)  is  a  slice  of 


.  ]5  - 


X  (a>i,  wj).  Re  {X  («!,  iDj)}  must  change  sign  somewhere  in  the  («lf  wj)  plane.  If  we 
assume  that  the  zeros  of  A  (z)  are  equally  likely  to  be  inside  the  unit  circle  as  outside, 
then  the  probability  of  A  (z  )  having  at  least  one  zero  outside  the  unit  circle  is 

P  «  HO Sf^'1  (14) 

since  a[n]  is  N -ft 2  points  long.  Since  this  condition  is  sufficient  but  not  necessary  for 

Re  {X  («j,  wj)}  to  change  sign,  p  represents  a  lower  bound  on  the  probability  that  a 

random  first-quadrant  sequence  would  satisfy  the  sign-change  constraint  of  Theorem  3. 

We  can  also  note  that  p  approaches  unity  very  rapidly  even  for  small  two-dimensional 

sequences.  For  example,  for  a  3x3  sequence,  p  -  09961,  and  for  a  4x4  sequence, 

p  =0999969.  For  a  64x64  sequence,  p  ~  1— 10“°°°. 

Thus,  we  have  shown  that  the  probability  of  a  random  first-quadrant  signal  satis¬ 
fying  the  constraints  of  Theorem  3  very  rapidly  approaches  unity  as  the  number  of 
points  in  the  signal  increases.  Although  this  does  not  guarantee  that  a  particular  signal 
in  some  particular  application  will  satisfy  the  constraints  of  Theorem  3,  it  does  show 
that  the  result  applies  to  a  broad  class  of  signals. 

22.  Unique  Specification  with  Samples  of  S,(uv  uj) 

Many  types  of  signals  or  functions  can  be  uniquely  specified  with  a  finite  set  of 
samples  evaluated  on  a  fixed  grid.  For  example,  the  DFT  is  a  method  of  sampling  the 
Fourier  transform  of  a  sequence  and  DFT  points  are  known  to  uniquely  specify  a 
finite  length  sequence.  Also,  in  phase-only  reconstruction  problems,  it  has  been  shown 
that  samples  of  the  phase  function  will,  for  the  most  part,  uniquely  specify  any 
sequence  which  is  uniquely  specified  by  its  complete  phase  function  [2,3,4]. 


Q 


In  the  problem  considered  here,  however,  the  information  in  Sx  (uj,  uj)  is  con¬ 
tained  in  the  exact  location  of  the  zero  crossings  and  this  information  is  lost  when 
Sx  (u»j,  u>j)  is  sampled.  From  another  point  of  view,  we  can  say  that  a  finite  set  of  sam¬ 
ples  of  S2  (<i>j,  (Oj)  contains  a  finite  number  of  bits  of  information  and  thus  cannot  be 
expected  to  uniquely  specify  a  signal  to  infinite  precision.  This  is  distinctly  different 
from  typical  sampling  problems  where  each  sample  is  of  (theoretically)  infinite  preci¬ 
sion  and  thus  does  not  contain  a  finite  number  of  bits  of  information.  (Note,  however, 
that  we  are  referring  strictly  to  theoretical  sampling  problems;  in  practical  applications 
signals  are  generally  represented  with  a  finite  number  of  bits,  and  it  may  be  possible 
for  a  signal  to  be  represented  to  sufficient  accuracy  with  a  finite  set  of  samples  of 
Sx  («,,  «^).) 

In  this  section,  we  shall  take  a  somewhat  different  approach  to  sampling 
Sx  («,,  <1*2).  Instead  of  using  the  values  of  £x(«*j,  <02)  over  a  fixed  grid,  we  shall  use  the 
location  of  points  on  the  zero  crossing  contours  of  Re  {, X  (<*)lP  012)}  as  the  'samples*. 
Since  these  contours  consist  of  an  infinite  number  of  points,  it  is  worthwhile  to  deter¬ 
mine  if  a  finite  set  of  such  points  will  uniquely  specify  a  signal.  Since  Theorem  1 
specifies  the  number  of  points  where  two  two-dimensional  polynomials  can  both  be 
zero,  we  can  use  this  theorem  to  establish  that  a  particular  number  of  arbitrarily- 
chosen  zero  crossing  points  is  guaranteed  to  be  sufficient  for  unique  specification.  We 
shall  also  show  that  this  number  of  points  may  not  be  necessary  for  unique 
specification;  in  particular,  if  the  zero  crossing  points  are  not  chosen  arbitrarily  but  are 
chosen  in  some  particular  way,  a  smaller  set  of  zero  crossing  points  can  be  sufficient 
for  unique  specification. 


.•  V  v  .• 
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Since  rectangular  regions  of  support  are  common  in  applications,  we  will  develop 
a  result  which  states  requirements  in  terms  of  such  regions.  The  result  could  be  easily 
modified  for  different  regions  of  support  or  could  be  applied  directly  to  a  problem 
involving  a  different  region  of  support  by  simply  assuming  a  rectangular  region  of  sup* 
port  large  enough  to  enclose  the  actual  region  of  support.  If  reference  to  a  region  of 
support  R  (N )  specifies  only  that  a  sequence  is  zero  outside  the  region 
-ft  s  nun2  ^ N ,  then  we  can  state: 


Theorem  6.  Let  x  [nx,  aj  and  y  aj  be  two-dimensional  sequences 
with  region  of  support  over  a  non  symmetric  half -plane  contained  in  R  (V  ). 
If  and  Yt  (*  i,  z  2)  are  nonfactorable  and 

Re  {X  (o>i»  02)}  -  Re  {Y  («,,  wj)}  =*  0  at  more  than  IfiJV2  distinct  points,  then 
x  l" t*  nJ  “  cy  [ji i,  nj  for  some  real  constant  c . 


Proof:  Recall  that  the  proof  of  Theorem  2  requires  stating  that  two  polynomials 


2j  Zj  Xt(z  1,^2)  and  zf  zj  Ye(z1,z2)  are  equal  to  within  a  scale  factor  given  that  they 
are  both  zero  at  an  infinite  number  of  points.  In  the  case  of  Theorem  6,  we  know  that 
zx  z$  X,(zlfZ2)  -zi  Yt(zx,z£  *0  at  more  than  lfitf2  points.  These  polynomials 
are  of  degree  4 N  and  thus,  by  Theorem  1,  can  have  at  most  lfitf 2  common  zeros.  Thus, 
z\  z2  Xt (z j, z 2)  *  czx  z*  Yg(zuz£  and  the  theorem  follows. 

As  we  did  with  Theorem  2,  we  can  modify  this  result  slightly  so  that  only  one 
sequence  is  required  to  have  a  nonfactorable  z-transform: 


Theorem  7.  Let  x  [/1  j,  aj  and  y  [alt  aj  be  two-dimensional  sequences 
with  region  of  support  over  a  nonsymmetric  half-plane  contained  in  R  ( N ). 
If  Xt(zx,zj)  is  nonfactorable,  x[V,AfJ#0,  and 
Re  {X  (<dj,  <1*2)}  =*  Re  {K(«i,  -  0  at  more  than  161V2  distinct  points,  then 

z  [a  1,  e2]  =  cy  [n j,  ej]  for  some  real  constant  c . 


"  »  ?  *•  m*  O  • 


The  proof  is  analogous  to  the  proof  of  Theorem  3. 

Note  that  for  this  result,  it  is  necessary  to  assume  the  nonfactorability  constraint 
of  Theorem  2.  If  a  sequence  satisfies  the  constraints  of  Theorem  S  but  not  the  con¬ 
straints  of  Theorem  2,  then  a  finite  set  of  zero  crossings  is  not  guaranteed  to  uniquely 
specify  the  sequence  since  these  zero  crossings  may  ail  correspond  to  the  same  factor. 
However,  a  similar  result  is  easily  developed  by  constraining  the  set  of  zero  crossings 
to  include  Nt  zero  crossings  of  the  /'*  factor,  where  each  Nt  is  chosen  in  accordance 
with  degree  of  the  factor. 

Next,  note  that  the  number  of  arbitrarily-chosen  zero-crossing  points  sufficient  to 
uniquely  specify  a  signal  is  somewhat  greater  than  the  number  of  points  in  the  signal. 
Specifically,  if  a  sequence  x  [nlf  nj  has  support  over  the  largest  possible  nonsymmetric 
half-plane  contained  in  R  (N),  then  it  has  2 N2  +2N  + 1  distinct  points.  According  to 
Theorem  6,x[n1,nJ\  is  uniquely  specified  with  m  zero  crossing  points  if  m  >  lfiAT2.  If 
x[it],n2]  is  real,  then  due  to  symmetries  in  the  Fourier  transform,  x  [»],  is  uniquely 
specified  with  m  >8N 2  points  if  these  points  are  chosen  over  an  appropriate  range  of 
frequencies.  This  means  a  sequence  consisting  of  p  real  points  is  uniquely  specified 
with  approximately  4 p  zero  crossing  frequencies.  If  a  sequence  has  support  over  only 
one  quadrant,  then  it  has  approximately  N2  points,  and  approximately  8 N2  zero  cross¬ 
ing  points  are  sufficient  for  unique  specification. 

Although  Theorem  6  states  that  a  particular  number  of  zero  crossing  points  is 
sufficient  for  unique  specification,  it  does  not  state  that  this  number  of  points  is  neces¬ 
sary  for  unique  specification.  In  particular,  as  we  show  next,  it  is  possible  to  specify  a 
signal  consisting  of  p  samples  with  p  -1  zero  crossing  points  if  the  zero  crossing  points 
are  not  chosen  arbitrarily  but  are  specifically  chosen  so  that  they  uniquely  specify  the 
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signal.  To  establish  this  result,  note  that  we  can  write  a  set  of  linear  equations  of  the 


form: 


COS  =  0 


where  each  equation  uses  a  different  pair  of  frequencies  («j,  <02)  on  a  zero  crossing 


contour,  ie,  where  the  equality  is  known  to  hold.  If  we  assume  that  x  [/tj,  nj  satisfies 


the  constraints  of  Theorem  7,  then  x  {V  ,N ]  =*  0,  so  we  can  substitute  x  [V  ,N ]  *  1  and 


obtain  a  non-zero  solution.  Thus,  if  x  [nlf  nj  consists  of  p  points,  these  equations  con¬ 


tain  p  -1  unknowns.  Although  we  have  not  shown  that  p  -1  equations  of  this  form  are 


guaranteed  to  have  a  unique  solution,  we  know  from  Theorem  7  that  if  a  sufficient 


number  of  equations  is  used  (say,  m  equations)  then  these  equations  are  guaranteed  to 


have  a  unique  solution.  If  we  have  m  equations  in  p  -1  unknowns  and  m  >  p  -1,  then 


some  of  the  equations  must  be  dependent  and  can  be  eliminated.  Thus,  it  is  possible 


(theoretically)  to  find  p  -1  independent  equations  from  the  set  of  m  equations,  and 


thus  the  corresponding  p  -1  zero  crossing  frequencies  are  sufficient  to  uniquely  specify 


x  [jij,  a  2].  This  result,  however,  does  not  suggest  a  practical  algorithm  for  choosing  the 


p  -1  zero  crossing  points  so  that  these  points  uniquely  specify  the  signal. 


3.  Application  to  Reconstruction  from  Zero  Crossings 


As  mentioned  in  the  introduction,  one  promising  application  for  the  results 


presented  in  section  2  involves  reconstructing  a  two-dimensional  signal  from  its  zero 


crossings,  or  more  generally,  from  its  threshold  crossings.  Considerable  research  has 


been  devoted  to  the  problem  of  reconstructing  one-dimensional  signals  from  zero 


crossings  (see  Requicha  [13]  for  a  survey.)  Most  of  this  work  has  concentrated  on  iden¬ 


tifying  types  of  one-dimensional  signals  which  have  a  sufficient  number  of  zero  cross- 
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ings  for  unique  specification.  The  most  recent  work  in  this  area  is  by  Logan  [14],  who 
states  a  set  of  conditions  under  which  one-dimensionai  bandpass  signals  are  uniquely 
specified  by  zero  crossings.  However,  no  corresponding  work  has  been  reported  on  the 
problem  of  reconstructing  two-dimensional  signals  from  zero  crossings.  In  addition, 
the  constraints  imposed  in  [14]  are  distinctly  different  from  the  constraints  to  be 
imposed  in  this  section. 

By  interchanging  the  roles  of  the  signal  and  transform  domains  in  Theorem  2, 
that  is,  by  exploiting  the  duality  of  the  Fourier  transform  in  a  straight-forward  way,  we 
can  develop  a  new  result  on  the  unique  specification  of  two-dimensional  signals  with 
zero  crossings.  Specifically,  if  a  continuous-time  signal  corresponds  to  the  Fourier 
transform  of  a  finite-length  discrete-time  sequence,  and  if  this  finite-length  sequence 
satisfies  the  conditions  of  any  of  the  theorems  developed  earlier,  then  the  signal  is 
uniquely  specified  to  within  a  scale  factor  by  its  zero  crossings.  To  illustrate  this  in 
more  detail,  consider  a  real,  band-limited,  continuous-time,  periodic  signal  /  (x ,  y  )  with 
periods  Tx  and  T2  in  the  x-  and  y-  directions,  respectively.  With  the  substitution  of 

variables  o>1  =  ^r~,  oij  »  -^L,  /  will  be  periodic  with  period  2ir  in  each  of  the  vari- 
M  T2 

ables  a)],  o> j  and  will  correspond  to  the  Fourier  transform  of  some  complex  conjugate- 
symmetric  finite-length  discrete-time  sequence  F  (nlf  nj): 
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expansion  of  f(x,y).  By  the  results  developed  earlier,  F(nlt /ij)  will  be  uniquely 
specified  by  sign  {f  (x  ,y )}  if  the  z-transform  of  F  (nu  is  nonfactorable.  The  z- 
transform  of  F  (n j,  nj)  can  be  obtained  from  the  right-hand  side  of  (16)  with  the  substi¬ 
tution  of  (z  f 1 ,  z  2_1 )  f°r  (c  Jm*),  (or  equivalently,  from  (17)  with  the  substitution 
-£ 'jr- 

of  (z  f 1 ,  z  2_1 )  for  (e  1  ,e  2 )  )  and  thus  represents  a  complex  extension  of  the  sig¬ 
nal  f(x,y),  where  f(x,y)  is  considered  to  be  known  over  the  unit  surface  in  the 
(z  j,  z  2)  complex  space.  Therefore,  the  z-transform  of  F(nun£  will  be  factorable  if 
and  only  if  /  (x ,  y ),  expressed  as  a  polynomial  in  equation  (17),  is  factorable. 

Let  us  state  our  dual  result  as  a  theorem: 


Theorem  8.  Let  f  (x  ,y)  and  g(x  ,y)  be  real,  two-dimensional, 
doubly-periodic,  continuous,  band-limited  functions  with 
sign  f  (x  ,y)  —  sign  g(x  ,y),  where  f(x,y)  takes  on  both  positive  and  nega¬ 
tive  values.  If  f  (x ,y)  and  g(x ,y)  are  nonfactorable  when  expressed  as 
polynomials  in  the  form  (17),  then  /  (x ,  y  )  ■  eg  (x ,  y  ). 


Proof:  This  theorem  follows  directly  from  Corollary  1  to  Theorem  2.  Since 
/ (x  ,y)  and  g(x,y)  have  Fourier  Series  coefficients  which  are  complex  conjugate  sym¬ 
metric,  have  finite  support,  and  have  nonfactorable  z-transforms,  they  satisfy  the  con¬ 
straints  of  Corollary  L  Thus,  the  Fourier  series  coefficients  must  be  equal,  and 
f(*>y)  =cg(x,y). 


Since  many  two-dimensional  signals  encountered  in  practice  are  not  periodic  but 
have  finite  support,  we  will  next  modify  Theorem  8  so  that  it  applies  to  these  signals. 
We  will  consider  the  case  where  /  (r ,  y  )  is  a  finite  segment  of  a  periodic  signal  satisfy¬ 
ing  the  constraints  of  Theorem  8.  For  example,  if  /  (x  ,y )  represents  one  period  of  a 
band-limited  periodic  function: 

/(*.y)=n/(x  +"2r2)  (jg) 

*1  *1 

then  it  is  possible  to  recover  /(x,y)  from  its  zero  crossings  provided  that  f(x,y) 
satisfies  the  constraints  of  Theorem  8,  even  though  f(x,y)  itself  is  not  band-limited. 
In  general,  it  is  not  necessary  for  the  duration  of  /  (x  ,y )  to  be  equal  to  one  period  of 
the  corresponding  periodic  function.  Thus,  /  (x  ,y  )  can  represent  a  finite  segment  of  a 
variety  of  different  periodic  functions.  In  order  for  f(x,y)  to  be  uniquely  specified 
by  its  zero  crossings,  we  only  need  one  periodic  function  containing  /(x,y)  to  be 
band-limited.  Specifically,  let  us  state  (see  Appendix  1  for  proof): 

Theorem  9.  Let  /(x,y)  and  g(x,y)  be  twodimensional  continuous 
functions  defined  over  the  same  known  region  R  of  finite  extent  with 
sign  f  (x  ,y)  *■  sign  g (x , y  ),  where  /  (x , y  )  takes  on  both  positive  and  nega¬ 
tive  values.  If /(x,y)  «/,(x,y)  in  R  and  g(x,y) -g,(x,y)  in  R  for  any 
periodic,  continuous,  band-limited  functions  /,(x,y)  and  gp(x,y)  which 
are  nonfactorable  when  expressed  as  polynomials  in  the  form  (17),  then 
/  (x  ,y)  =  eg (x , y )  for  some  positive  constant  c . 

Note,  however,  that  Theorems  8  and  9  do  not  imply  that  an  arbitrary  (non¬ 
periodic)  band-limited  continuous  function,  or  a  finite  segment  of  one,  is  uniquely 
specified  by  its  zero  crossings.  These  results  would  be  equivalent  to  developing  a  result 
similar  to  Theorem  2  for  continuous-time  signals. 


.23. 


r 

l 


r 


l* 


Although  Theorems  8  and  9  were  developed  from  Corollary  1  to  Theorem  2,  the 
other  extensions  developed  for  Theorem  2  can  also  be  useful  here.  Of  particular 
interest  to  us  is  the  fact  that  the  definition  of  sign  information  can  be  extended  to 
include  crossings  of  an  arbitrary  threshold  instead  of  just  zero  crossings.  This  result  is 
important  since  in  a  number  of  applications,  particularly  image  processing,  signals  are 
constrained  to  be  nonnegative  and  thus  have  no  zero  crossings.  Nevertheless,  it  might 
be  desirable  to  recover  such  a  signal  from  knowledge  of  the  points  where  it  crosses 
some  particular  level.  The  same  procedure  used  to  establish  Theorem  4  can  also  be 
used  to  show  that  Theorems  8  and  9  can  be  extended  to  include  reconstruction  from 
sig n(f  ( x ,  y  )-p)  provided  that  3  is  chosen  so  that  the  signal  actually  crosses  the  thres* 
hold  p.  In  this  case,  if  we  have  sign(f  (x ,  y  )-fJ)  =  sign(g  (x ,  y  )-p),  then 
f(x,y )  -P=c(g(x,y)  -  P)  for  some  c . 

4.  Reconstruction  Algorithms 

Having  established  that  particular  classes  of  signals  are  uniquely  specified  by  some 
partial  information,  it  is  of  interest  to  develop  algorithms  for  recovering  the  original 
signal  from  this  information.  One  common  approach  to  developing  algorithms  for 
reconstruction  from  various  forms  of  partial  information  is  to  develop  an  iterative 
algorithm  which  alternately  impose  constraints  in  the  space  and  frequency  domains. 
Another  approach  is  to  express  the  solution  as  a  set  of  simultaneous  linear  equations. 
In  this  section,  we  will  discuss  each  of  these  methods  and  present  experimental  results 
obtained  with  each  method. 


a 
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41  Iterative  Alforithm 


The  class  of  iterative  algorithms  mentioned  above  can  be  applied  to  the  problem 
of  reconstruction  from  Fourier,  sign  information  by  imposing  the  correct  sign  of  the 
real  part  of  the  Fourier  transform  in  the  frequency  domain  and  imposing  the  known 
region  of  support  in  the  space  domain.  It  can  also  be  applied  to  the  problem  of  recon¬ 
struction  from  sign  if  (x  ,y)}  by  imposing  the  correct  sign  (perhaps  with  respect  to 
some  threshold)  in  the  space  domain  and  the  band-limited  constraint  in  the  frequency 
domain.  For  the  sake  of  clarity,  in  this  section  and  the  one  that  follows  we  will  pri¬ 
marily  refer  to  the  reconstruction  from  Ss  (<alt  uj);  however,  it  should  be  realized  that 
the  same  algorithms  can  be  applied  to  reconstruction  from  sign  {/  (x ,  y )}  by  inter¬ 
changing  the  roles  of  the  two  domains. 

Since  knowledge  of  the  exact  points  of  discontinuity  is  necessary  for  the  signal  to 
be  uniquely  specified  (in  either  problem),  the  convergence  of  an  iterative  algorithm  to 
the  correct  solution  necessarily  depends  upon  the  use  of  the  exact  zero  crossing  points. 
Thus,  an  algorithm  for  reconstruction  from  5,  (wj,  wj  which  uses  a  DFT  and  thus 
makes  use  of  only  those  values  of  5,  (ult  012)  corresponding  to  DFT  points  cannot  be 
guaranteed  to  converge  to  the  correct  solution.  However,  it  can  be  shown  that  the 
continuous-frequency  version  of  the  algorithm  (that  is,  a  similar  algorithm  imposing 
the  correct  value  of  Sx  (wj,  <1*2)  for  all  frequencies  and  thus  using  an  actual  Fourier 
transform  and  not  a  DFT)  will  converge  to  the  correct  sequence.  It  is  also  possible  to 
show  that  the  sampled-frequency  version  of  the  algorithm  will  converge  to  a  sequence 
which  satisfies  both  the  space  and  frequency  domain  constraints  (provided  such  a  solu¬ 
tion  exists),  although  the  solution  is  not  unique.  These  resuits  can  be  developed  within 
the  theory,  of  projections  onto  convex  sets,  as  was  used  in  [15]  to  establish  the 
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convergence  of  a  number  of  different  signal  reconstruction  algorithms.  Specifically, 
the  results  developed  in  [15]  apply  directly  to  this  problem  provided  the  constraints  in 
each  domain  are  imposed  in  such  a  way  as  to  be  projections  onto  convex  sets;  the 
details  will  be  presented  in  Appendix  2. 

Once  the  theoretical  properties  of  a  reconstruction  algorithm  have  been  deter¬ 
mined,  it  is  important  to  empirically  determine  the  effectiveness  of  the  algorithm  in 
recovering  an  actual  signal.  In  particular,  it  is  worthwhile  to  determine  if  a  practical 
size  DFT  limits  the  set  of  solutions  to  a  sufficiently  small  set  and  if  convergence  (or  a 
good  approximation)  can  be  obtained  with  a  practical  number  of  iterations.  It  is  also 
worthwhile  to  investigate  the  effect  of  using  different  initial  estimates  in  the  iteration. 
Experimentally,  we  have  found  that  if  the  DFT  size  used  is  at  least  4  times  the  signal 
size  and  an  initial  estimate  is  chosen  which  in  some  sense  resembles  the  original  signal, 
then  reasonable  results  can  be  obtained  with  a  small  number  of  iterations.  Specifically, 
we  have  had  good  success  when  using  an  initial  estimate  formed  from  the  correct  one 
bit  of  Fourier  transform  phase  and  a  Fourier  transform  magnitude  which  is  the  aver¬ 
age  of  a  number  of  unrelated  images.  Although  a  large  number  of  iterations  are 
required  for  the  algorithm  to  converge  to  a  sequence  satisfying  the  time  and  frequency 
domain  constraints,  the  improvement  of  image  quality  after  the  first  20  iterations  or  so 
is  somewhat  negligible  even  if  the  image  at  this  stage  does  not  satisfy  the  frequency 
domain  constraints  at  every  point.  An  example  is  included  in  Figure  1,  where  we  show 
the  original  image  (a)  and  the  image  reconstructed  from  one  bit  of  Fourier  transform 
phase  (b).  In  this  example,  the  original  image  is  64x64  points,  256x256  DFTs  were 
used,  and  the  results  shown  were  obtained  with  25  iterations. 


42.  Linear  Equation  Method 


As  mentioned  earlier,  it  is  also  possible  to  express  the  solution  to  the  problem  of 
reconstruction  from  one  bit  of  phase  as  a  set  of  linear  equations.  Thus,  another  possi¬ 
ble  reconstruction  algorithm  would  involve  solving  the  set  of  equations: 

22  x  [nlf  n  J  cos  +<1)202)  s0 

•1  «i  '  ' 

where  each  equation  uses  a  different  pair  of  frequencies  («j,  wj)  for  which  the  equality 
is  known  to  hold. 

As  mentioned  earlier,  if  we  assume  that  x  [nj,  nj  satisfies  the  constraints  of 
Theorem  7,  then  xflV.JV]  #0  and  we  can  substitute  x {V  ,N ]  *  1  and  obtain  a  non-zero 
solution.  If  the  number  of  equations  is  chosen  in  accordance  with  this  corollary,  then 
these  equations  are  guaranteed  to  have  a  unique  solution.  However,  this  number  is 
significantly  higher  than  the  number  of  equations  usually  required  in  practice.  Our 
experience  is  that  for  a  sequence  with  p  nonzero  points,  p  -1  equations  in  p  -1  unk¬ 
nowns  will  generally  have  a  unique  solution,  although  the  results  are  particularly  sensi¬ 
tive  to  numerical  errors  in  the  values  of  the  zero  crossing  frequencies.  If  the  values  of 
and  a>2  are  obtained  to  four-  or  five-digit  accuracy  and  if  the  number  of  equations 
used  is  only  slightly  greater  than  the  number  of  unknowns  and  a  least-squares  solution 
is  obtained,  results  indistinguishable  from  the  original  signal  can  be  achieved.  We  have 
used  this  procedure  to  reconstruct  a  number  of  different  images  of  varying  sizes  which 
satisfy  the  constraints  of  the  results  developed  earlier  and  have  always  successfully 
recovered  the  original  image  provided  enough  equations  were  used. 

An  example  of  results  obtained  with  this  method  is  shown  in  Figure  2,  which 
shows  the  original  image  (a)  and  the  image  reconstructed  by  solving  the  above  equa- 


tions  (b).  In  this  example,  the  original  image  is  15  (120  points)  and  122  equations  in 

119  unknowns  were  used.  The  real  parts  of  the  Fourier  transforms  of  these  images  are 
shown  in  Figures  3  (a)  and  (b).  If  the  image  of  Figure  3  (a)  is  considered  to  be  the  ori¬ 
ginal  image,  then  Figure  3  (b)  represents  an  image  reconstructed  from  its  zero  cross¬ 
ings.  For  comparison,  an  image  showing  only  the  zero  crossings  of  the  image  in  Figure 
3  (a)  (in.  Figure  3  (a)  quantized  to  one  bit)  has  been  included  as  Figure  3  (c). 

5.  Conclusions 

In  this  paper,  we  have  developed  conditions  under  which  two-dimensional  signals 
can  be  reconstructed  from  Fourier  transform  sign  information.  We  have  also  shown 
that  these  conditions  apply  to  a  broad  class  of  signals  by  showing  that  if  the 
coefficients  of  the  signal  are  random,  the  probability  of  a  sequence  satisfying  these 
conditions  very  rapidly  approaches  unity  as  the  number  of  points  in  the  signal 
approaches  infinity.  From  the  basic  result,  we  have  developed  a  number  of  extensions 
to  different  types  of  signals,  different  constraints,  and  different  types  of  sign  informa¬ 
tion.  We  have  also  applied  these  results  to  the  problem  of  reconstructing  a  two- 
dimensional  signal  from  the  zero  crossings  of  the  signal  itself.  In  addition,  we  have 
discussed  some  possible  algorithms  for  reconstruction  from  one  bit  of  Fourier 
transform  phase  or  from  the  signal  zero  crossings.  Examples  of  images  reconstructed 
with  these  algorithms  have  also  been  included. 


Appendix  L  Proofs 


This  appendix  contains  proofs  of  some  of  results  presented  earlier. 

Proof  of  Corollary  1  to  Theorem  2.  This  proof  follows  the  proof  of  Theorem  2, 
with  a  few  modifications.  First  of  all,  since  the  sequences  are  permitted  to  be  complex, 
the  real  part  of  the  Fourier  transform  corresponds  to  the  conjugate  symmetric  com¬ 
ponent  of  a  sequence,  ian  (keeping  the  notation  x«[*i,  "J) 


With  this  substitution,  the  proof  of  Theorem  2  is  then  directly  applicable  to  the  case  of 
causal  complex  signals,  with  the  exception  that  Im  {x  [OjOj}  is  not  recoverable  from 
x,[n1,ji2]-  For  conjugate  symmetric  signals,  we  note  that  x,[/Ji, nj  nji  and  thus 

Theorem  2  and  its  proof  are  directly  applicable  with  the  substitution  of  xjnj,  nj  for 

» 

*,[«i,  "21- 

Proof  of  Theorem  3.  Following  the  proof  of  Theorem  2,  we  know  that  since 
Sz toj)  s  S, («!,  aij),  X0(zvz£  and  Yt(zltz£  must  contain  a  common  factor.  Since 


Xc(z  zj)  is  nonfactorable,  then 


if  *€>[!«!, n J, 


Yt (z i, z 2)  ~ X t(z i, Z2)  F (z\,z 2)  for  some  real  symmetric  factor  F(zj,z 2).  Equivalently, 

where  f[nv nj  is  even.  Thus,  *0  for  some 

values  nl>  N  or  n2>  N ,  violating  the  constraints  of  the  theorem.  Thus, 
*  [« 1* #* 2]  "21- 

Note  that  since  the  convolution  involves  xf  (n  j,  nj  and  not  z  [n,,  nj,  knowledge 
that  x  [0, 0]  #  0  is  not  sufficient  to  guarantee  uniqueness,  but  knowledge  that 
x[/V  J  *  0  is  sufficient.  The  information  required  is  the  exact  size  of  the  region  of 
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support  of  xt[nit  nj*  which  is  the  same  regardless  of  the  value  of  x  [0, 0J. 

Proof  of  Theorem  4.  Note  that  with  the  definitions  given  in  the  theorem  state¬ 
ment, 

*(«1,«2)  "Re  {X  fa, i^)eJa}  (A2) 

Y  («!,  «j)  -  Re  {Y  («!,  u£eJa}  -  0 

Then,  we  note  that  £/*(<»},  U]) -5i  fa,  w^,  and  since  «2>, 

5_-  (oilf  taij)  =  («j,  wj).  Since  i[«i,  nj  and  >[«i, /ij  are  complex  conjugate  symmetric 

sequences  and  satisfy  the  constraints  of  Corollary  1  to  Theorem  2,i[n1,n2]  =  cy  [n^nj 
for  some  positive  constant  c.  Then,  with  some  algebra,  xfnj,  nj  =  cy[/t1,  itj  for 
0*1,02)  *  (0,0),  and  x[0,0]  cos  a  -  0  **c(y  [0,0]  cos  a  -  0)  for  some  positive  constant 
c. 

Proof  of  Theorem  5.  Recall  from  the  proof  of  Theorem  2  that  if  Xt(zx>z^)  and 
K#(zx,z  2)  have  common  zero  contours  on  the  unit  surface  then  they  must  have  a  com¬ 
mon  factor.  We  will  assume  that  x  [ji^  hJ  *  cy  [a  j,  02]  and  attempt  to  reach  a  contrad¬ 
iction.  For  convenience,  let  us  assume  that  there  is  some  irreducible  factor  of 
X,(z1,z2)  which  is  not  a  factor  of  Yt(zx,zj).  First  of  all,  note  that  if  this  factor, 
denoted  F{?x,zj),  is  complex,  then  F*(zx,z£  will  also  be  a  factor  of  Xt(zx,zj)  and 
thus  X,(zx,z 2)  will  contain  a  real  symmetric  factor  F(zl,z£F*(zl,z£  which  is  nonne¬ 
gative  everywhere,  violating  the  constraints  of  the  theorem.  If  this  factor  F(zx,z£  is 
real  but  nonsymmetric,  then  F(zf1,z2”1)  will  also  be  a  factor  since  X,(zt,Z2)  «  sym¬ 
metric.  In  this  case,  Xt(zx,z£  will  contain  a  real  symmetric  factor 


,  v  v  >■  "*.  '  k  .  •  —  w  h\.  *  .  *  -  «  w  r  •  /  .•  /  .•  •  .* 
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F  (z  j,  zj)  must  be  real  and  symmetric,  and  since  according  to  the  theorem  hypothesis,  it 
has  both  positive  and  negative  values  and  has  multiplicity  one,  then  we  must  have 
Sx  (<!>],  <i>2 )  *Sy  (at},  u)j)  for  some  values  of  (wj,  ooj),  and  we  have  reached  a  contradic¬ 
tion.  Thus,  there  cannot  be  any  factor  of  JTff(z1,22)  which  is  not  a  factor  of  Yt(zuz£ 
and  thus,  x  [n  j,  aj  *  h  i»  "J* 

Proof  of  Theorem  9z  First  we  note  that  if  fp(x,y)  and  gp(x  ,y)  are  not  simply 
periodic  replications  of  /  (x ,  y  )  and  g  (x ,  y  ),  then  fp(x,y)  and  gp  (x ,  y  )  may  contain 
zero  crossings  that  cannot  be  obtained  from  the  zero  crossings  of  f(x,y)  and  g(x ,y). 
Thus,  Theorem  9  does  not  quite  follow  directly  from  Theorem  8.  However,  as  we  saw 
in  the  proof  of  Theorem  6,  unique  specification  in  terms  of  zero  crossings  does  not 
require  knowledge  of  all  the  zero  crossing  points;  it  requires  only  a  finite  set  of  them. 
In  this  case,  as  long  as  f  (x  ,y)  contains  both  positive  and  negative  values,  it  will  con- 
tain  at  least  one  zero  crossing  contour  with  an  infinite  number  of  points.  Thus 
fp(x,y)  -gp(x  ,y)  =  0  at  an  infinite  number  of  points,  and  using  arguments  taken 
from  the  proofs  of  Theorems  2  and  0,fp(x,y)  mcgp(x,y).  Since  f  (x  ,y)  and  g(x,y) 
are  both  known  to  be  defined  over  the  same  region  R  and  over  R ,/(x,y)  =fp{x,y) 
and  g(x,y)  =  gp(x  ,y),then  /  (x  ,y)  -  eg(x,y ). 


Appendix  2.  Convergence  of  Iterative  Reconstruction  Algorithm 

In  this  appendix,  we  establish  the  convergence  of  our  algorithm  for  reconstruction 
from  one  bit  of  Fourier  transform  phase  by  using  the  theory  of  projection  onto  convex 
sets,  using  the  approach  presented  in  [15]  for  establishing  the  convergence  of  a  variety 
of  signal  reconstruction  algorithms.  Specifically,  the  result  we  shall  be  using  is  as  fol¬ 
lows*: 

Theorem  AL  [1546]  Let  H  be  a  Hilbert  space,  G  be  a  composition  of 
projection  operators  onto  closed  convex  sets,  at  least  one  of  which  is  finite¬ 
dimensional,  and  G*  denote  the  intersection  of  these  sets.  If  G*  is 
nonempty,  then  for  all  x  €  H ,  the  sequence  Gmx  converges  to  a  point  in  G* . 

The  results  developed  in  [15]  include  showing  that  a  wide  variety  of  constraints 
often  imposed  in  signal  reconstruction  algorithms  can  be  imposed  in  such  a  way  that 
the  transformations  will  satisfy  the  constraints  of  Theorem  L  In  fact,  the  constraints 
used  in  our  algorithm  for  reconstruction  from  Ss  (tuj,  uj)  differ  from  some  of  those  dis¬ 
cussed  in  [15]  in  only  trivial  ways.  The  basic  approach  shall  be  repeated  here  although 
the  mathematical  details  shall  be  omitted. 

To  show  that  Theorem  1  applies  to  our  iterative  algorithm,  we  must  first  carefully 
define  the  transformations  applied  in  the  time  (or  space)  and  frequency  domains  at 
each  iteration  so  that  they  can  be  characterized  as  projection  operators  onto  closed 
convex  sets.  First  let  us  note  that  although  the  continuous-frequency  variables  (<■»},  uj) 
will  be  used  throughout  the  discussion  below,  the  properties  of  the  transformations  in 
the  time  and  frequency  domains  apply  equally  well  if  the  discrete  frequencies 

corresponding  to  a  DFT  are  used.  The  only  difference  is  that  in  the  continuous- 
Theorem  A1  ia  a  weak  fora  of  the  results  in  [1346]  hut  it  is  sufficient  for  our  purposes  here. 


frequency  case,  the  set  G*  will  contain  exactly  one  sequence  (provided  the  proper  con¬ 
straints  are  satisfied),  whereas  in  the  discrete-frequency  case,  the  set  G*  will  contain  an 
infinite  number  of  sequences. 

Let  h  denote  the  sequence  we  are  trying  to  recover  which  is  known  to 

satisfy  the  constraints  of  Corollary  1  to  Theorem  2  and  is  thus  uniquely  specified  by 
Sk  (oij,  ojj).  Let  T*  denote  the  set  of  sequences  which  satisfy  the  time  domain  con¬ 
straints: 

“0  for  ^  P/f]  (A3) 

xV/,N]=h\?f,N] 

and  F*  denote  the  set  of  sequences  which  satisfy  the  frequency  domain  constraints: 

Re  (X  («,,  wj)}  a  0  if  Sk  (»,,  wj)  ■  1 

Re  (X  («!,  «2)}  s  0  if  u»2)  -  ~1  (  ' 

Note  that  the  set  T*  is  finite-dimensional  even  if  the  space  H  includes  infinite- 
length  signals.  Also  note  that  the  definition  of  the  set  F*  is  not  precisely  the  same  as 
stating  SI  («,,  wj)  -  Sh  («,,  wj);  the  difference  occurs  if  Re  {H  («,,  ovj)}  <  0  and 
Re  (X  (wj,  «j)}  =  0  at  some  point.  It  is  necessary  to  use  the  definition  of  F*  given  in 
(4)  in  order  for  F*  to  be  a  closed  set  (a  set  which  includes  its  limit  points.) 

Next  we  will  define  the  operators  7  and  F  to  be  projections  onto  the  sets  T*  and 
F * .  For  7  and  F  to  be  projections,  we  need: 

|  |7x  -*  1 1  s  |  \y  -x  1 1  for  all  y  €  7*  (A5) 

and 

1 1  Fx  -x  J  |  s  1 1  y  -x  | )  for  all  y  €  F*  (Afi) 

Thus  the  operator  7  we  need  to  impose  time  domain  constraints  is: 
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r# 


0 


7  *2)] 


0  (*1,1*2)  *(N,N) 

h(N  ,N)  (1*1, 1*2)  *  (N  ,JV) 

0  »i,»2  > Af 


(A7) 


or  in  words,  simply  substituting  the  known  values  of  h[n1,n2].  The  operator  F  we 
need  to  impose  frequency  domain  constraints  is: 


F l* 0*„ *2)]  ~*/pr  («i.  <•*)],  (A8) 

where 


F/pTK.wj)] 


*  (•»!,  «2)  ^  -S,  («,,  «2)  -  5*  («1,  W2) 

Im  {X  («!,  ^2)}  otherwise 


or  in  words,  keeping  the  imaginary  part  constant  and  setting  the  real  part  to  zero  if  its 


sign  is  incorrect. 

Next,  we  express  our  iterative  algorithm  in  the  form 


xkH~Gxk  (A9) 

where  G  =TF  is  a  composition  of  projection  operators.  Then,  by  Theorem  1,  the 

sequence  in  equation  (9)  will  converge  to  a  point  in  G*,  that  is,  a  sequence  which 
satisfies  the  time  and  frequency  domain  constraints.  Thus,  if  satisfies  the  con¬ 

straints  of  Theorem  3,  and  the  iteration  imposes  the  correct  Sh  (a»j,  02)  for  all  frequen¬ 
cies  (in.,  actual  Fourier  transforms  are  used),  then  G*  contains  exactly  one  sequence, 
and  the  iteration  must  converge  to  that  sequence.  If  Sk  («,,  <o 2)  is  sampled  (in,  a  DFT 
is  used),  then  the  iteration  must  converge  to  a  sequence  in  G* ,  in,  a  sequence  which 
satisfies  the  time  and  frequency  domain  constraints,  although  this  solution  is  not 
unique  and  the  solution  actually  obtained  depends  on  the  initial  estimate. 


[» 
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Figure  2.  Reconstruction  with  Linear  Equation  Method 


(a)  original  image 


(b)  recovered  image 


Figure  3.  Reconstruction  from  Zero  Crossings 
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